X 



Asymptotic Inverse Problem 
for Almost-Periodically Perturbed 
Quantum Harmonic Oscillator 



Q ■ Alexis Pokrovski 

O 
(N 

00 



Abstract 



Let {^„}^Q be the spectrum of —-^+x +q{x) in L (M), where q is an even almost- 
periodic complex-valued function with bounded primitive and derivative. Suppose that 



the asymptotic approximation to the first asymptotic correction A/i„ = — //^ -|- 
o{n~^), /x^ = 2n -|- 1 is known. We prove the formula that gives the frequencies and 
the Fourier coefficients of q in terms of A/i„. 



\Q [ AMS Mathematics Subject Classifications (1991): 34L20 (Asymptotic distribution 

O ■ of eigenvalues, asymptotic theory of eigenfunctions) ; 

■ 81Q15 (Perturbation theories for operators and differential equations). 
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a 

> ; 1 Introduction and main result 



Consider the operator describing perturbed quantum harmonic oscillator 

A = -— + x^ + q{x) in L\R) (1.1) 

with the perturbation g(a;) from the class i3 = {q : ||^'||oo + ||Q||cxd < oo}, where Q(a:) = Jf^qdi 
and II • lloo denotes the norm in L°°(R). It was proved in jlj that the spectrum {/i„}J!^Q of A 
has the asymptotics /i„ = /i° + /i^ + 0(n~3), where = 2n + 1 and //^ = 0{n~i). 

For the perturbations that are sum of almost-periodic and decaying terms we study the 
problem of recovering of the almost-periodic part from the first asymptotic correction /x^. 
Specifically, we consider the perturbations 

1 

q = p + rGB, p E ,p(—x) = p{x) and ||r||Bi = lim — / \r(x)\dx = 0, (1.2) 

T^oo 2T J _x 
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where is the Besikovitch space of almost-periodic functions [5] (closure of trigonometric 



polynomials Ylk=o^kG^'^''^ , real, in the norm || ■ HbO- It is sufficient to recover p in terms 
of its Fourier transform [2]. Here is the main result. 

Theorem 1.1. Let {A/i„}^^ approximates the first asymptotic correction to the spectrum 
of the operator M.l]) . n~^) : 

A/i„ = /i„ - /i° + o(?2"3). (1.3) 

Then the spectrum and the Fourier coefficients of the almost-periodic part p can be re- 
covered from the relation 

1 1 f'^ 

lim — AfinGuixn,XL){xn+i - Xn) = hm — / p(t)cosutdt, u ^ 0, (1.4) 
^ T Jo 

where x„ = a/T^ = v^^rT+T, Gy{x,T) = —x f^-^^^^-, fu,T(t) = r]{t — T)cosz/t and 
rj G C^(]R) is a smoothed step function such that r]{t) = 1 for x G (— c>o, —1], rjit) = for 
X G [0, oo) and r]'{0) = 0. 

Asymptotic inverse spectral problem for quantum harmonic oscillator with slowly de- 
caying perturbation was considered by Gurarie 0. He studied the operator p.lj) with real 
q{x) ~ o,m cosUmX for |a;| —>■ oo, where a > and the sum in the numerator is finite. 

The approach in [2| is based on the spectral asymptotics 



^l/4+a/2 



+ 0{—=), g(a;) = const — ^ cos(ti;ma; — 7r/4) 



which exhibits linear relation between the leading asymptotic terms of q and [in- However, 
the technique of [S] does not cover the case a = 0. 

We consider just this case in a slightly more general setting (almost-periodic functions 
vs. finite trigonometric sums). Our method also allows complex- valued q. Technically, the 
result is based on the recent proof [T] of the spectral asymptotics 

/^n = /x° + 7^/ g(Vp^sin^9)ci^9 + 0(n-i), for g G i3. (1.5) 

Thus the proof of Theorem 11.11 follows from the asymptotic behavior of the integral in p.5|) . 
which is analyzed in Lemmas 12.11 and 12.31 

2 Properties of the Schlomilch integral 

The integral in the spectral asymptotics ()1.5|) is the Schlomilch integral jl] 

9,{x) = - r^\4xsm^)d^=- r q4^) = iqix) + qi~x))/2 (2.6) 

Jo Jo VX^ - t^ 



evaluated at the points Xn = + 1. In the next Lemma we estimate the integral and 
its derivatives. Then in Lemma 12.31 we prove similar estimates for the inverse Schlomilch 
integral. (We could not find in the literature the results of these Lemmas for the specific 
class Using the two Lemmas, we prove Theorem ll.il 
Everywhere below C denotes an absolute constant. 

Lemma 2.1. Let f G B and g{x) = J^^'^ f{xsm'd)d'd. Then 

\g^^)\^C^13^^±M^, \g'(x)\^C^^^^^iM^, x>0, (2.7) 
\/l + X Vl + X 

where F{x) = f dt. 

Proof. For a; ^ 1 the result is evident, so we consider only the case a; > 1. Using the 

■ f{xt)di 



change of variables t = sin??, we write g{x) = tttt^ split it as 



f{xt)dt ^ f{xt)dt 



where e = 1/x. Using the notation (. . .)[ for ^(. • •) and choosing the primitive F of /, 
satisfying F(x(l — e)) = 0, we have 

l-e X / / ~ , , t=l-e l-e =, 



A = ^ / (H^t) j . 7^ = + / iSfFi ) ■ Therefore, 



|/,|<2^(^l + /^j<cM^(l + -L). (2.9) 

Now we substitute the estimate {hi ^ ||/||oo /q'' :J = 2||/||oo/yi and (Q into (EH). This 
gives the first inequality in (|2.7p . We prove the second one in a similar way, writing 

' j'^j' J' r'tnxt)dt r tf'[xt)dt 

We integrate by parts in I[, choosing the primitive f{xt) = f{xt) — /(x(l — e)). This gives 

f{xt) dt 
(l-t2)3/2' 



r - -i r^^" ^^""^^ '^^ hence 
-'i ~ X h (i_t2)3/2, neuce 



X 



^C*^^ ( 1 + ^ ) . (2.11) 



We substitute the estimate I/2I ^ C\\f'\\oo/y/x and ()2.11|) in ()2.1()|1 . This gives the second 
inequality in ()2.7j) . ■ 

Remark 2.2. The rate of decay as x ^ 00 in p.7|l cannot be improved. (The example 
/(x) = cosx gives the Bessel function Jq.) 
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Lemma 2.3. Let T > 2, (p,ip" G L°°{[0,T]) and ^p{T) = 0. Then the equation ip{t) = 
f It^ J^^a^^^a has the unique solution g{x) = —x ^^^^'^ /^'^ ^ ^ [0' '^'"^^ ^^^^ 

|^7(a;)| ^C(||<^|U + ||v^'||oo)v^, for x > 1. (2.12) 
//, in addition, f'{T) = 0, then 

\g'{x)\ ^ Ciy'lU + ||v^"||oo)v^, for x > 1. (2.13) 
Proof. In terms of g{x) = and (p{t) = (p{\/t) the equation on g becomes the Abel 
or^nofi^Ti ,~('7-'i = 1 ^!/^ ■ solution for absolutely continuous (p is g{s) = ^^rz^ — 



'^J^l^ (^^^ Ch.l, §2 of ini). Using ip{T) = 0, we obtain the required formula for g. 
Consider ()2.12j) . For x G [T — 1,T] the inequality follows from the direct estimate 
|£{£)| ^ rj dt _ p 3, g rg 2- _ 1] ^rite 

«<^> A^rm. /.^rm (2.14) 



x^ 



and integrate I2 by parts. We have I2 



fit) 



t=x+l 



- L+i ^it)aiJW^ Therefore, 



|/2|^^^ + ||^||oo-i=iI="°° (2.15) 
Vl + 2x - t=x+i VI + 2x 

Now we substitute the estimate ^ ||v^^||oo f^^^ y/J-x'^ ^ ||v''||oo\/2/x and ()2.15|) into 
(jmH). This gives ^U^ . as required. 

Next consider ()2.13p . By g'{x) = x{g{x) / x)' + g{x) / x and ()2.12|) . it is sufficient to estimate 
x{g{x) / x)' . Using ^p'{T) = 0, we obtain 

For X e [T - 1,T] we have \x{g{x)/xy\ ^ W'{t)\dt ^ 2^/2y\\^^. For x G 

[0, T - 1] write 

+ 4 /;^r^, 4^r^ (2.17, 

and take the integral for by parts. We have J2 = "■^^=^^=4 + Ix+i^'(^)§t^/W^ 
Hence, using ^yp=f ^ we obtain 

(l + a;)||^'|U , uju r d i-t) 



I/2I ^ ^ + ll'^'ll- / ^^^=^^^ ^ 2y\UVT+^. (2.18) 

V 1 + 2a; Jx+i 0** vt^ - x^ 

Now we substitute the estimate |/[| ^ WWooJ^'^^ ^ llf^'llooVl + 2x and (ITT^ into 

(|TT7|1 . This gives (EH- ■ 
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Remark 2.4. Note that the condition (p{T) = is necessary for ()2.12|) . If v^(T) 7^ 0, 
then g{x) is unbounded due to the non- integral term in the inversion formula for the Abel 
equation. The term is 0(^==) for x ] T. Similarly, the condition ^'{T) = is necessary for 



(j2.13|) . If v^'(T) 7^ 0, then g'{x) is unbounded as a; t T due to the non-integral term ^ 

in (inni). 

Proof of Theorem II. IL Compare with the asymptotis ()1.5j) . It is clear that 



^fJ'n = gq{xn) + o{n i) , x„ = a/ /i^ = \^2n + 1, (2.19) 

where the Schlomilch integral Qq is given by ()2.6|) . The proof is based on the fact that the 
set {xn}^=Q becomes arbitrarily dense as n — > 00, so that Riemann sums 

Sx„+i<r Gu{xn, T)gq{xn){xn+i-Xn) approximate ^ Gy{x, T)gg{x)dx, provided the integ- 
rand is smooth enough. Using the inversion formula for the Schlomilch integral, we choose 
Gi, such that the last expression tends to ^ cos utdt as T ^ 00. 

By Lemma f2. 31 we have 



G,ix,T)gq{x)dx = ^^/§=f-) = u{t)^u,T{t) dt. (2.20) 

Let us show that the integrand in the left-hand side of ()2.20|) is sufficiently smooth. By 
Lemma 12.11 

|^,(x)|^C^^±^, \g'g{x)\^GM^^±M^, x^O, (2.21) 

where C is an absolute constant. Similarly, since (pi^^T satisfies the hypothesis of Lemma [2.31 
uniformly in T ^ a; 

|G',(x,T)KC(l + z/)v^, |fa(a;,T)| ^C(l + z/)2y^, x ^ 1, (2.22) 

where we used ||v5i/,r||oo ^ G, Wf^xWoo ^ G{1 + z/)^. Therefore, for any fixed z/ the function 

hrix) '= G,y{x,T)gq{x) and its x-derivative are uniformly bounded for x ^ 1 and T ^ x. 
Hence, for T ^ cxd we have 

- hT{xn){xn+i- Xn) - - j hrix) dx = - 0{-) dx = 0{—^) ^ 0, (2.23) 

where we used Xn+i — Xn = 0{x^^). Now, by ()2.23|) . ()2.19p and the first estimate in ()2.22|) . 
lim — Afj,nGu{xn,T){xn+i - Xn) = lim — / Gu{x,T)gq{x)dx. (2.24) 

Next we divide ()2.20|) by T and take the limit T oo. By ()2.24j) and lim qj^ipjy^x dt = 
\im p{t) cos ut dt, this gives ()1.4|) .B 

Remark 2.5. We have (uxy^Guix, oo) = ^ = Joi^x), where Jq is the Bessel 
function (see e.g. 0). 
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